3D hybrid computations for streamer discharges and production of run-away electrons 
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We introduce a 3D hybrid model for streamer discharges that follows the dynamics of single 
electrons in the region with strong field enhancement at the streamer tip while approximating the 
many electrons in the streamer interior as densities. We explain the method and present first results 
for negative streamers in nitrogen. We focus on the high electron energies observed in the simulation. 
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Streamers are a fundamental mode of electrical break- 
down of ionizable matter when a strong voltage is ap- 
plied; they are the first stage in the evolution of sparks 
and lightning. Streamers are ionized plasma channels 
that grow into a non-ionized medium due to the self- 
enhancement of the electric field at their tips (see Fig. [2]). 
In this high field region, the electron energy distribution 
is very far from equilibrium and can have a long tail 
at high energies [H, [E HIi which makes the streamer a 
plausible candidate for the generation of so-called run- 
away electrons [E IE Hi- Such very energetic electrons 
subsequently can produce X-rays and 7-rays through 
Bremsstrahlung, therefore they may explain X-ray bursts 
and flashes observed during thunderstorms 0, [1] , rocket 
triggere d lig htning [9(| and spark development in the lab- 
oratory [la III Hill- 

Streamer dynamics is mostly modeled by a fluid (or 
density) model 0, EE EE El, EE El as this approxi- 
mation is computationally most efficient and able to de- 
scribe the main characteristics of the streamer discharge 
in a qualitative way. However, it obviously cannot trace 
the single particle dynamics. To follow the distribu- 
tion of positions and velocities of individual electrons 
- and therefore density fluctuations, run-away effects 
and excited molecular levels — , a particle (or Monte 
Carlo) model [E HE HH is required that follows individ- 
ual electrons and their elastic, inelastic and ionizing colli- 
sions with the background of abundant neutral molecules. 
However, the particle model is not suitable to study the 
single electron dynamics either, because the increasing 
number of electrons eventually renders computational 
power and storage unaffordable, while a super-particle 
approach causes numerical heating and stochastic arti- 
facts 

We therefore here introduce a hybrid streamer model 
for full three-dimensional calculations. It uses the natu- 
ral structure of the streamer: the particle model is ap- 
plied in the most dynamic and exotic region with rela- 
tively few electrons and high local electric field, i.e., in 
the ionization front, and the many slow electrons inside 
the streamer channel are left to the fluid model. How to 
implement the spatial coupling of density and fluid model 
in the one-dimensional case, was presented in [E Q and 
illustrated in Fig. 1 in both papers. Here the method is 



extended to 3D, and first results for a negative streamer 
in nitrogen at standard temperature and pressure are pre- 
sented. We first discuss model and numerical implemen- 
tation and then the physical results. 

Particle and extended fluid model and Poisson solver. 
The particle and the fluid model were compared quan- 
titatively in E|, Q for negative streamers in nitrogen for 
electric fields from -190 Td to ~750 Td, or from 50 
kV/cm to 200 kV/cm at standard temperature and pres- 
sure. An important finding in [l( was that the electron 
and ion density in fluid or particle model start to dif- 
fer when the field exceeds ~190 Td, and, more impor- 
tantly, that this relative difference largely increases with 
increasing electric field. In the reasons of this density 
discrepancy are discussed, and it is shown that the fluid 
model had to be extended by a gradient expansion to op- 
timally approximate the particle model. This model has 
the form 

dne w • c dn P c n\ 

W + v ' Je = 5 ' -dT = s > (1) 

j e = —fi(E)En e — D(E) • Vn e , (2) 
S = n(E)a(E){E n e + k 1 (E)E-Vn e ), (3) 

where n e and n p are electron and ion density, respec- 
tively, j e is the electron flux and S is the nonlocal source 
term with a density gradient expansion parameterized by 
ki, fi represents the mobility and D is the diffusion ten- 
sor, and E = (E x , E y , E z ) and E are the electric field and 
its strength. The electric field is calculated with a fast 
Poisson solver: the 3D fishpack subroutine [IE HI]. 

Differential cross-sections. The cross sections for the 
relevant collisions in the 3D particle model are taken from 
the SIGLO database [25j for incident electrons with ener- 
gies up to 1 keV. Above 1 keV, the Born approxima- 
tion [26J is used for elastic collisions, a fit formula in (2?| 
is implemented for the electronically exciting collisions 
and the Born-Bethe approximation [IE HI is used for 
ionizing collisions. The electron transport coefficients, 
reaction rates and the average energies are generated in 
particle swarm experiments [l], Q ; they agree well with 
the Boltzmann solver (BOLSIG+) [IE HE] when in both 
cases isotropic scattering and equal energy sharing in 
ionizing collisions is assumed. The scattering method 
derived by Okhrimovskyy et al. 31] is implemented for 
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elastic and exciting collisions. Opal's empirical fit [12] is 
implemented for the energy splitting in ionizing events, 
where incident electrons with high energies are likely to 
keep most of their energy. 

The spatial coupling of fluid and particle model — more 
precisely, the position of the model interface as a function 
of the maximal field and the construction of the buffer re- 
gion — was already discussed in @, 0] for planar fronts. 
When the 3D streamer is decomposed into many nar- 
row parallel columns oriented in the propagation direc- 
tion (as detailed further below), this coupling can be ap- 
plied in each of these columns. However, new problems 
arise due to the complexity of the 3D geometry: i) The 
model interface in 3D is never planar, but depending on 
the used criterion, it is either smoothly curved or even 
strongly fluctuating; a fluctuating model interface will 
create large buffer regions and dramatically increase the 
computational cost. Since in small grid cells, the electric 
field is smooth while the electron density can fluctuate 
heavily, the position of the model interface is determined 
here through the electric field rather than through the 
electron densities. More precisely, in the results shown 
below, the model interface in each column is placed where 
the field is E = 0.85 E + with E + being the maximal field 
ahead of the front within the column. This criterion en- 
sures that the relative error for the electron densities in 
the streamer stays below 3% for all fields E + [3]. The 
large region at the sides of the streamer that stays non- 
ionized, is treated by the particle model, ii) A direct 
contact of particle and fluid model without a buffer re- 
gion can cause electron leaking, and hence loss of mass 
and charge. Therefore the buffer region has to be con- 
structed carefully not only at the ionization front, but 
also in the lateral directions. Details on the model inter- 
face and the buffer region are given after introducing the 
structure of the simulation results. 

A hybrid streamer simulation. Figs. HHH show different 
aspects of the same simulation. It is a negative streamer 
in nitrogen at standard temperature and pressure. It 
propagates through a gap of 1.18 mm between two pla- 
nar electrodes; the applied voltage is 11.8 kV which cor- 
responds to a background field of 100 kV/cm or 372 Td. 
The simulation starts with 100 electrons and ions sit- 
ting 0.05 mm away from the cathode. They are ini- 
tially followed by the pure particle model, and the hy- 
brid model is introduced at time 0.32 ns when the num- 
ber of electrons reaches 1.5 x 10 7 in a manner discussed 
further below. The simulations are carried out on a uni- 
form grid of 256 x 256 x 512 grid points with the cell 
length Ax = Ay — Az = 2.3 /xm and with time step 
At = 0.3 ps, the numerical procedure for particle and 
fluid model are described in [2, H, [H . 

Fig. Q] shows the electron energy distribution at the 
moment when the simulation switches from pure particle 
to hybrid computations; the curved model interface is 
also marked. The figure shows that the region with high 
mean electron energies is covered by the particle model 
and the low energy part with many electrons is left for 




FIG. 1: Local mean energy of the electrons at time 0.32 ns 
when the simulation switches from pure particle to hybrid 
simulation; the buffer regions between particle and fluid 
regime are indicated with a dark curve line. 



the fluid model; here the fluid model is both efficient and 
appropriate. Furthermore, the particle model is applied 
in all regions of low to vanishing electron density ahead 
and at the sides of the streamer; here the particle model 
is both more correct and also more efficient than the fluid 
model. 

Fig.[2]shows the electric field E z in the direction of the 
background field at times 0.36 ns, 0.45 ns and 0.54 ns of 
the simulation (cf. [13, [33[ for a more extended discussion 
and more plots of the streamer evolution in fluid approx- 
imation). The location of the buffer region is marked 
in red. The fluid model is applied within the red lines 
and the particle model is applied the large outer region 
where the field enhancement region is always included. 
As E z is large, the buffer region in z-direction should be 
2 or even 3 cells long to obtain a stable electron flux at 
the model interface [3j ; this procedure is applied in each 
column where one column is one row of cells in the z di- 
rection. In the x- and y-direction, one cell is long enough 
for the buffer region since the radial electric field is much 
smaller, but to prevent electron leaking from the particle 
region directly to the fluid region, 2 cells are used. 

In practice, the hybrid simulation is very efficient in 
approximating the majority of the electrons by densities 
and in following the streamer much longer than the pure 
particle model. Specifically, at the times 0.36 ns, 0.45 ns 
and 0.54 ns shown in Fig. H only 12%, 7% and 3% of 
the electrons are followed individually. Nevertheless, the 
figure shows that in the region with the highest electric 
field the single electrons are followed. We remark that 
this simulation costs 43 hours on a normal desktop (Intel 
Quad2 CPU, 8 Gb RAM). 

Run-away electrons. The electron-nitrogen collision 
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FIG. 2: The electric field E z after time 0.36 ns, 0.45 ns and 0.54 ns of the simulation with model interface and buffer regions 
marked in red. All quantities are shown on two orthogonal planes that intersect with the 3D structure. The full system size is 
0.59 x 0.59 x 1.18 mm 3 ; only the dynamically interesting regions are shown. 
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FIG. 3: The electrons with energy above 200 eV at time steps 
0.36 ns, 0.45 ns and 0.54 ns. The maximal field E + at these 
times is 160, 220, and 290 kV/cm. 



frequency is maximal for electron energies of about 
200 eV, beyond that energy they have a chance to run 
away as the friction decreases when the energy increases 
further. Fig. [3] therefore shows only the electrons with 
energy above 200 eV at the same three time steps as 
Fig. O Electrons with e > 200 eV start to appear when 
the maximal field E + reaches 160 kV/cm. But these 
electrons lose their energy almost immediately again. As 
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FIG. 4: The maximal electric field strength E + , the number 
of electrons 7V e >2oo cv with energy larger than 200 eV and 
the highest electron energy e max as a function of time. The 
maximal electron energy exceeds 1 keV at t~0.54 ns. 



the maximal field E + increases further during streamer 
propagation, both the number and the energy of the high 
energy electrons increases. Although most electrons still 
very quickly lose their energy, a few ones are able to ac- 
celerate further, and at time 0.54 ns, electrons with en- 
ergy above 1 keV are observed. When the streamer later 
approaches the upper anode, the field increases further, 
also due to the proximity of the electrode, and electron 
energies up to 3.5 keV are seen. 

Fig. 2] analyzes the situation further. Plotted is the 
maximal electric field strength E + , the number of elec- 
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trons with energy above 200 eV, and the highest electron 
energy. Until approximately 0.2 ns, the maximal field 
equals the background field, i.e., the system is in the 
avalanche phase and no energetic electrons are present. 
After time 0.3 ns, the maximal field enhancement in- 
creases more than linearly in time, after 0.36 ns the first 
electrons above 200 eV appear, and after 0.45 ns their 
number and energies increase massively. Large fluctua- 
tions in the maximal electron energy as a function of time 
can be seen; there is not one electron that runs away, but 
many are being accelerated on average. Given the dis- 
tance of 0.25mm that the front crosses between times 
0.45 ns and 0.54 ns, the maximal electrostatic energy of 
the background field is » 2.5 keV; over this distance elec- 
trons accelerate from 0.2 to 1 keV. 

During the time interval from 0.36 ns to 0.54 ns, the 
field at streamer head is enhanced to 1.5 to 3 times of 
the background field. These fields can accelerate elec- 
trons beyond the maximum of the electron-neutral fric- 
tion force (cf. Fig. 2 in [5] or Fig. 9 in [6]) of 200 eV. 
Electrons in the energy range of several hundred eV get 
well ahead of the front, but many of them do not fully 
run away. As they get from the region of enhanced elec- 
tric field to the region where the field decays ahead of the 
front while inelastic and ionizing scattering is still con- 
siderable, they are trapped and create many new small 



avalanches ahead of the ionization front [34], [35, 36] . The 
electrons with energies of several keV are likely to keep 
accelerating even in the lower background field ahead of 
the streamer p| [33] , but in the present simulation they 
rapidly reach the anode and disappear. 

We have presented a 3D hybrid model for streamers 
that reliably can follow the single electron dynamics in 
the high field region of the streamer head at moderate 
computational costs, and that can observe electrons be- 
ing accelerated to over 1 keV. Electrons with energies 
above 200 eV appear when the field enhancement at the 
streamer head exceeds 160 kV/cm or 600 Td. The en- 
ergetic electrons can run out of the streamer head and 
relax somewhat ahead of the ionization front creating 
new avalanches; in this way they can create local front 
jumps and increase the mean velocity of the front. The 
investigation of streamers in air rather than in nitrogen 
will be subject of future studies, as well as the question 
whether streamers powered by higher voltages, e.g., in 
the corona of lightning leaders, can accelerate electrons 
into the relativistic range of MeV energies. 
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